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Abstract 

The polarized antiquark distributions in the proton can be measured by studying spin 
asymmetries in vector boson production in longitudinally polarized proton-proton collisions. 
The STAR and PHENIX experiments at BNL RHIC have reported first observations of single 
spin asymmetries in VF^-production most recently. We compute the QCD corrections to 
single and double spin asymmetries, taking account of the leptonic decay of the boson 
and of restrictions on the kinematical acceptance of the detectors. The QCD corrections 
have only a small impact on the asymmetries, such that a reliable extraction of the polarized 
antiquark distributions can be envisaged once more precise measurements are made. 
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1 Introduction 



Owing to a large variety of measurements, the unpolarized parton distributions are known to 
high precision [1-5]. Determinations of the polarized parton distributions [6,7] have to rely 
on a much smaller data set, consisting mainly of inclusive and semi-inclusive polarized deep 
inelastic scattering measurements on fixed targets, augmented with single inclusive hadron and 
jet measurements from RHIC. As a result, the polarized quark distributions in the proton are 
determined accurately, while constraints on the polarized antiquark and gluon distributions are 
much less stringent. Semi-inclusive polarized deep inelastic scattering provides some information 
on the antiquark distributions of different flavors, although this extraction is not free from 
ambiguities, since it requires the quark-to-hadron fragmentation functions as input. 

A very clean measurement of the polarized antiquark distributions in the proton can be 
made from spin asymmetries in electroweak vector boson production in polarized proton-proton 
collisions [8]. Owing to the parity-violating coupling of the and bosons, non- vanishing 
spin asymmetries are already obtained in the collision of one polarized and one unpolarized 
proton (single spin asymmetries), in addition to the more common double spin asymmetries. 
Measurements of these asymmetries are possible at the BNL RHIC collider, and both the STAR 
and PHENIX experiments have recently reported a first observation of single spin asymmetries 
in production [9,10]. At present, these data are not yet sufficiently precise (and do not cover 
a large enough kinematical range) to provide competitive constraints on the polarized antiquark 
distributions. However, future polarized proton runs at RHIC will improve upon this situation 
and are likely to yield precision measurements of the spin asymmetries in production. 

To include these data in global fits of polarized parton distributions at next-to-leading order 
(NLO), the NLO QCD corrections must be taken into account in the evaluation of the asymme- 
tries. These corrections are known already for a long time for the fully inclusive asymmetries 
in gauge boson production [11, 12] and for the asymmetries differential in the gauge boson ra- 
pidity [13,14]. Since VF^-boson production is observed only through the decay into a lepton 
and a neutrino, the VF^-boson rapidity cannot be measured directly. By accurately determining 
the missing transverse momentum due to the unobserved neutrino and measuring the lepton 
transverse momentum and rapidity, the 1^^-boson rapidity can in principle be inferred from 
kinematical constraints. With the limited coverage of the STAR and PHENIX detectors, such 
a determination is however not possible. Instead, one measures the asymmetries as a function 
of the kinematics of the observed final-state lepton: lepton rapidity yi and lepton transverse 
momentum pT,i- These correlate closely with the VF^-boson rapidity, but the interpretation of 
the resulting asymmetries in terms of the polarized parton distributions is less straightforward. 

It is the aim of this paper to derive the NLO QCD corrections to the single and double 
spin asymmetries in production as function of the observed lepton rapidity and transverse 
momentum. For the single spin asymmetries, these corrections were recently obtained by de 
Florian and Vogelsang in [15], and our calculation provides an independent validation of their 
results. Higher-order QCD effects in the asymmetries, including resummation of large loga- 
rithmic corrections, were studied previously in [16] and are implemented in the widely used 
RHICBOS program. 
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This paper is structured as follows: in Section 2, we introduce the single and double spin 
asymmetries in VF^-boson production and discuss their interpretation in terms of polarized 
parton distributions. Section 3 briefly reviews the technical aspects of the NLO calculation, 
while Section 4 discusses the numerical impact of the NLO corrections. We compare the NLO 
results to the recent STAR and PHENIX measurements in Section 5 and conclude with Section 6. 

2 Single and double spin asymmetries in vector boson produc- 
tion 

The single and double spin asymmetries in vector boson production processes are constructed 
from the helicity-dependent production cross sections 

where hi^2 = (+) — ) denote the helicities of the two incoming hadrons. The coordinate frame is 
defined such that hadron 1 moves in positive rapidity direction. 

The unpolarized cross section is the average over all initial state helicities: 

da = ^ (d(7++ + d(7+- + du'+ + da—) , (1) 
while the singly (only hadron 1 polarized) and doubly polarized cross sections are given by: 

dA(7L = ^ (df7++ + d(7+- - d(7-+ - da") , (2) 
dAaLL = \ (da++ - da+- - du'^ + da~) . (3) 

Prom these, single and double spin asymmetries are constructed as functions of kinematical 

variables of the observed final-state particles. Previous studies on NLO corrections to spin 
asymmetries in vector boson production [14] focused on the rapidity of the vector boson y, 
which is itself not directly measurable in production, due to the unobserved neutrino in the 
final state. 

For the partonic interpretation of the asymmetries, the vector boson rapidity is very instruc- 
tive, since at leading order (LO) it is directly related to the momentum fractions of the partons 
probed in both hadrons: 

Restricting to only up and down quarks, one thus finds very simple leading-order expressions 
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for the asymmetries: 



-l^u{x\)d{xl) + /\d{xl)u{xl) 
u{x\)d{xl) + d{x\)u{xl) 
^d{x\)u{x%) + Au(xg)d 
d{x^{)u{x%) + u{x^^)d{xl 
^u{x\)M{x%) + A(i(x°)AK(x°) 
u(xO)(J(xO) + (i(x?)M(.xO) 



^ ^ Ad(x°)An(x^) + A^(xO)Ad(xO) 

At RHIC with \fS = 500 GeV, these asymmetries are thus sensitive on xi^2 ~ 0.05. At large and 
positive rapidity y, {x^ > X2), the single spin asymmetries Al are dominated by the first term in 
the numerator, since Ag(x5, Q"^) ^ ^q{x\, Q^) for large x\. The second term in the numerator 
is dominant for large and negative y, {x\ < x%), since q{x2-iQ^) 3> ([{x^^Q"^) for large x%. Given 
that for the x-range probed by RHIC, the polarized quark distributions are substantially larger 
than the polarized antiquark distributions, and their sum is well constrained from polarized 
inclusive deep inelastic scattering [6], the double spin asymmetries can be used for a reliable 
extraction of An(x) and A(J(x). 

With only the decay lepton being observable, the lepton rapidity yi and transverse momentum 
Pt^i are the more appropriate kinematical variables. Measurements of the transverse mass, as 
commonly carried out in vector boson production at the Tevatron, are not feasible at RHIC due 
to the non-hermetic detector coverage, which prevents a reliable reconstruction of the missing 
transverse momentum carried by the neutrino. 

We shall thus focus on the single and double spin asymmetries 

dAfjL/dy; dAo-Li/dyj 

and 

dAo-L/dpT,; . , . _ dAo-LL/dpr,; 

Polarized and unpolarized cross sections in these asymmetries can be expanded in QCD pertur- 
bation theory in the strong coupling constant a^. In these cross sections, the exchanged 
boson is not required to be on-shell, its finite width is taken into account in its propagator. 



3 Next-to-leading order corrections 

The next-to-leading order expressions for the asymmetries are obtained by including the 0{as) 
corrections in numerator and denominator. The computation of these corrections is a standard 
task in higher order calculations, carried out in dimensional regular ization in d = 4 — 2e dimen- 
sions. It requires to combine the virtual one-loop corrections to the process qq' — > I'^u 
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with the real radiation contributions qq' W^g l^i'g and qg W^q' l^uq'. The real 
radiation contributions develop infrared singularities if the final-state gluon becomes soft, or 
if the final-state parton becomes collinear with one of the incoming partons. These infrared 
real radiation singularities can be extracted analytically. They cancel in the final expression 
for the cross section when combined with the infrared-divergent virtual contributions and with 
mass factorization counterterms for the incoming parton distributions. Various techniques for 
the treatment of infrared singular real radiation exist at NLO, based either on introducing 
subtraction terms [17-19] or on a slicing of the final-state phase space [20,21]. 

In our calculation, we use the phase space slicing method [20, 21] for the real radiation con- 
tributions. In this method, the real radiation phase space is split into resolved and unresolved 
regions by introducing a cut-off parameter s^i^ on all Mandclstam invariants. The resolved re- 
gions are integrated numerically, while the unresolved regions are integrated analytically, giving 
rise to the infrared divergent terms. The phase space slicing contributions from the unresolved 
regions are universal and depend only on the type of partons involved in the unresolved config- 
uration [20]. These are computed first in a hypothetical kincmatical situation with all partons 
in the final state, and then continued to the true kinematical situation by using a crossing func- 
tion [21] for each initial state parton. The polarized crossing functions are obtained from the 
unpolarizcd expressions in [21] by substituting the d-dimensional unpolarized one-loop splitting 
functions by their polarized counterparts [22,23]. 

Our calculation was validated by comparing the NLO unpolarized differential cross sections 
da against MCFM [24] and the singly polarized differential cross sections dAcr^ against CHE [15] 
(using the same set of parameters and setting renormalization and factorization scales = 
for the total cross sections and fj? = (M^^r + p|,^)/4 for the distributions, as in [15]). Full 
agreement is found for both quantities, thereby providing a crucial independent cross-check 
of [15]. 

4 Numerical results 

Polarized proton-proton collisions are studied at RHIC at different centre-of-mass energies. For 
vector boson production, only the high-energy runs at s/S = 500 GeV are of relevance, such that 
we use this setup for our studies. For the unpolarized parton distribution functions, we use the 
MRST2002 NLO set [25]; polarized parton distributions are taken from DSSV [6]. The DSSV 
set is the result of a global fit to spin asymmetries in inclusive and semi-inclusive deep inelastic 
scattering as well as single inclusive hadron and jet production at RHIC. In its extraction, the 
unpolarized cross sections entering the asymmetries are computed from the MRST2002 NLO 
set, such that the DSSV spin-dependent parton distributions should be used consistently with 
this set. Likewise, the value of the strong coupling constant should be consistent with what is 
used in the parton distributions, such that we take as{Mz) = 0.1195 [25]. Renormalization and 
factorization scales are taken to be equal, and are fixed to /ir = = Mw- The NLO scale 
dependence of the single and double spin asymmetries ^^(y) and ALL{y) was studied in detail 
in [14]. It turns out to be very small, largely due to cancelations between the polarized and 
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Figure 1: Next-to- leading order differential cross sections for single inclusive lepton production 
at RHIC, ^/S = 500 GeV with pT,i > 20 GeV and \yi\ < 2. Solid: unpolarized cross section, 
dashed: singly polarized, dotted: doubly polarized. 

unpolarized cross sections in the asymmetry. 

The mass and width are taken [26] as Mw = 80.399 GeV and Fw = 2.085 GeV, and the 
Fermi coupling as Gp = 1.16637 • 10"'^ GeV^^. By using the Fermi coupling constant as input 
parameter, electroweak corrections to the cross sections are minimized. Neglecting incoming 
bottom quarks, we use the CKM matrix elements \ Vud\ = |Ks| = 0.975 and |V^s| = \ Vcd\ = 0.222. 

As default for our numerical studies, we use the following cuts on the final-state lepton: 
Pt,i > 20 GeV and \rii\ < 2. The unpolarized and polarized NLO cross sections with these cuts 
are displayed in Figure 1 as function of either yi or px^i- From these results, it can be anticipated 




that the single spin asymmetries arc substantially larger than the double spin asymmetries. The 
transverse momentum distribution of the lepton shows the characteristic peak at in the 

unpolarized and polarized cases. To study the behavior of the polarized cross sections in more 
detail, the single and double spin asymmetries Al and All are more suitable than the cross 
sections themselves. 

Figures 2 and 3 show these asymmetries as function of yi and pT,i- In these figures, the 
LO asymmetry is obtained by evaluating polarized and unpolarized cross sections to LO, while 
both cross sections are evaluated to NLO for the NLO asymmetry. Since this figure aims to 
quantify the effect of the parton-level NLO corrections, both LO and NLO asymmetry are 
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Figure 3: Single and double spin asymmetries as function of pT,i at NLO (solid) and LO (dashed). 



evaluated with the same NLO sets of parton distributions. We observe that Aid/i) and All^jji) 
are only very mildly affected by the NLO corrections, and that the corrections do not change 
the shape of these asymmetries. The corrections to Al{pt,i) and All{pt,i) are also very small 
and uniform for px^i < Mw/2, while they become more pronounced for p^i > M\v/2. This 
behavior can be easily understood from kinematical considerations. At LO, the boson is 
produced at vanishing transverse momentum, such that leptons above pT,i = M\y/2 must come 
from off-shell production. The real radiation contributions at NLO lead to final states with 
a finite transverse momentum of the W^, which can thus be near its mass shell. With this new 
parton- level contribution appearing at NLO, the NLO unpolarized and polarized cross sections 
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Figure 4: Individual partonic contributions to the single and double spin asymmetries as function 
of ui- Unpolarizcd cross section is always evaluated to NLO, polarized cross section evaluated 
to NLO (solid), LO (dashed), NLO without qg subprocess (dotted). 



for px^i > are much larger than their LO counterparts, and the asymmetry is modified 

substantially. 

Compared to the results of [15], we observe slight discrepancies in ^L(y;) for large \yi\. These 
are entirely due to the different choice of renormalization and factorization scales: /v,^ = in 
our evaluation, compared to fj? = (M^ + j»|,|)/4 in [15]. With the same scale choice (and the 
same electroweak parameters), we find full agreement with [15]. 

The numerical importance of the different parton-level contributions to the single and double 
spin asymmetries is displayed in Figures 2 and 3. All asymmetries in these figures are normal- 
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Figure 5: Individual partonic contributions to the single and double spin asymmetries as function 
of pt^i- Unpolarized cross section is always evaluated to NLO, polarized cross section evaluated 
to NLO (solid), LO (dashed), NLO without qg subprocess (dotted). 



ized to the unpolarized cross sections evaluated at NLO. Comparing the results obtained with 
Act evaluated only at LO with the LO asymmetries in Figures 2 and 3, one concludes that the 
smallness of the NLO corrections to the asymmetries is due to a cancelation of the NLO correc- 
tions to unpolarized and polarized cross sections in the asymmetry. The NLO corrections to the 
polarized cross sections are moreover dominated by the qq subprocess. Consequently, a future 
extraction of polarized quark distributions from the Aiiyi) and ALL{yi) will not be affected by 
the uncertainty on the polarized gluon distribution. 
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Cuts 




A T ( PXD 1 


At ihO) 


At fNLOl 


STAR 


25 GeV < pT,i < 50 GeV; 

\yi\ < 1 


w+ 


-0.27 ±0.10 ±0.02 


-0.347 


-0.348 


w- 


0.14 ±0.19 ±0.02 


0.127 


0.123 


PHENIX 


30 GeV < pT,i < 50 GeV; 
\yi\ < 0.35 


w+ 


-0.861°:?° 


-0.342 


-0.342 


w- 


0.881°:^? 


0.107 


0.103 



Table 1: Single spin asymmetries measured in production at RHIC by STAR [9] and 
PHENIX [10] compared to LO and NLO QCD predictions. 



5 Comparison with RHIC data 

The single spin asymmetries Al in H^^-boson production have been measured recently at RHIC 
by the STAR [9] and PHENIX [10] experiments. Both measurements are based on detecting 
only electrons and positrons at large transverse momentum. Cuts are imposed on the lepton 
transverse momentum and rapidity, and the measurement is performed in a single bin. Especially 
the PHENIX measurement relies on a very limited kinematical coverage, such that and 
boson production cannot be disentangled, since the second lepton from the decay is often 
outside the coverage. The Z°-boson contributions to the single spin asymmetries at NLO were 
investigated in [15]. 

In Table 1, we summarize the experimental cuts and compare the experimentally measured 
asymmetries with the theoretical predictions at LO and NLO. We observe that the predictions 
for the asymmetries are perturbatively very stable, with practically no shift observed for the 
production, and only a modification of four per cent for the W~ production. The asymmetries 
observed by STAR arc in good agreement with the theoretical prediction, while PHENIX typ- 
ically obtains asymmetries with a larger magnitude than expected, although with large errors. 
Both measurements have clearly demonstrated the existence of single spin asymmetries in 
production, but not yet attained sufficient precision to provide meaningful constraints on the 
polarized antiquark distributions. 

Future more precise data on the single spin asymmetries, as well as first data on the double 
spin asymmetries will provide very important input to the determination of the polarized parton 
distributions. 

6 Conclusions 

Spin asymmetries in gauge boson production in polarized proton-proton collisions are currently 
measured at RHIC [9,10]. These asymmetries provide a direct measurement of the polarized 
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light antiquark distributions, which were only loosely constrained from semi-inclusive polarized 
deep inelastic scattering up to now [6,7]. In this paper, we derived the NLO QCD corrections 
to single and double spin asymmetries in boson production in polarized proton-proton 
collisions as function of the measurable final-state lepton variables. Our results for the single 
spin asymmetries confirm a recent calculation of de Florian and Vogelsang [15]. 

We show that the NLO corrections have only a very moderate impact on the asymmetries, 
and that the corrections are largely independent on the at present poorly constrained polar- 
ized gluon distribution. Their inclusion is nevertheless mandatory in future global NLO fits of 
polarized parton distributions for the sake of internal perturbative consistency. 
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